Mechanism for the Intermittent Route to Strange Nonchaotic Attractors 
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Intermittent strange nonchaotic attractors (SNAs) appear typically in quasiperiodically forced 
period-doubling systems. As a representative model, we consider the quasiperiodically forced logis- 
tic map and investigate the mechanism for the intermittent route to SNAs using rational approxi- 
mations to the quasiperiodic forcing. It is thus found that a smooth torus is transformed into an 
intermittent SNA via a phase-dependent saddle-node bifurcation when it collides with a new type 
of "ring-shaped" unstable set. Besides this intermittent transition, other transitions such as the 
interior, boundary, and band-merging crises may also occur through collision with the ring-shaped 
unstable sets. Hence the ring-shaped unstable sets play a central role for such dynamical transitions. 
Furthermore, these kinds of dynamical transitions seem to be "universal," in the sense that they 
occur typically in a large class of quasiperiodically forced period-doubling systems. 

PACS numbers: 05.45.Ac, 05.45.Df, 05.45.Pq 



Recently much attention has been paid to the study of 
quasiperiodically forced systems because of the generic 
appearance of strange nonchaotic attractors (SNAs) [Q. 
SNAs were first described by Grebogi et al. || and since 
then have been extensively investigated both numerically 
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imcntally |L7). SNAs exhibit some properties of regu- 
lar as well as chaoic attractors. Like regular attractors, 
their dynamics is nonchaotic; like typical chaotic attrac- 
tors, they exhibit fractal phase space structure. Further- 
more, SNAs are related to Anderson localization in the 
Schrodinger equation with a spatially quasiperiodic po- 
tential [fl8f , and they may have a practical application in 
secure communication fL9| . Therefore, dynamical tran- 
sitions in quasiperiodically forced systems have become 
a topic of considerable current interest. However, their 
mechanisms are still much less understood than those in 
unforced or periodically forced systems. 

Here we are interested in the dynamical transition to 
SNAs accompanied by intermittent behavior, as reported 
in [[l2|. As a parameter passes a threshold value, a 
smooth torus is abruptly transformed into an intermit- 
tent SNA. Near the transition point, the intermittent dy- 
namics on the SNA was characterized in terms of the av- 
erage interburst time and the Lyapunov exponent. This 
route to an intermittent SNA is quite general and has 
been observed in a number of quasiperiodically forced 
period-doubling maps and flows (e.g., see |l3|, Q). It 
has been suggested that the observed intermittent 
behavior results through interaction with an unstable or- 
bit. However, in the previous work, such an unstable 
orbit was not located, and thus the bifurcation mecha- 
nism for the intermittent transition still remains unclear. 

In this paper, we study the underlying mechanism of 
the intermittency in the quasiperiodically forced logis- 
tic map M H which is a representative model for the 



quasiperiodically forced period-doubling systems: 

_ f x n+ i = (a + e cos2Tr9 n )x n (l - x n ), , , 
' 1 Qn+i =On+u (mod 1), 1 > 

where a; S [0, 1], 6 G S , a is the nonlinearity parameter 
of the logistic map, and oj and e represent the frequency 
and amplitude of the quasiperiodic forcing, respectively. 
We set the frequency to be the reciprocal of the golden 
mean, uj — (V5 — l)/2- The intermittent transition 
is then investigated using the rational approximations 
(RAs) to the quasiperiodic forcing. For the case of the 
inverse golden mean, its rational approximants are given 
by the ratios of the Fibonacci numbers, LOk = Fk—i/Fk, 
where the sequence of {Fk} satisfies Fk+i = Fk + Fk-i 
with Fq = and Fi = 1. Instead of the quasiperiodi- 
cally forced system, we study an infinite sequence of pe- 
riodically forced systems with rational driving frequen- 
cies LOk- We suppose that the properties of the original 
system M may be obtained by taking the quasiperiodic 
limit k — > oo. Using this technique we observe a new 
type of invariant unstable set, which will be referred to 
as the "ring-shaped" unstable set in accordance with its 
geometry. When a smooth torus (corresponding to an 
ordinary quasiperiodic attractor) collides with this ring- 
shaped unstable set, a transition to an intermittent SNA 
is found to occur. 

We also note that the quasiperiodically forced logis- 
tic map M is noninvertible, because its Jacobian de- 
terminant becomes zero along the critical curve, Lq — 
{x = 0.5, 8 E [0,1)}. Critical curves of rank k, 
(k = 1,2,.. .), are then given by the images of L , [i.e., 
Lfe = M k (L )]. Segments of these critical curves can be 
used to define a bounded trapping region of the phase 
space, called an "absorbing area," inside which, upon 
entering, trajectories are henceforth confined pOfl . It is 
found that the newly-born intermittent SNA fills the ab- 
sorbing area. Hence the global structure of the SNA is 
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FIG. 1: (a) Phase Diagram in the a — e plane. Regular, 
chaotic, SNA, and divergence regimes are shown in light gray, 
black, gray (or dark gray), and white, respectively. To show 
the region of existence (gray) of the intermittent SNA occur- 
ring between T (light gray) and the chaotic attractor region 
(black), a small box near (a, e) = [3.38, e*(= 0.584 726 781)] is 
magnified. Through interaction with the ring-shaped unsta- 
ble set born when passing the dashed line, typical dynamical 
transitions such as the intermittency (route a) and the inte- 
rior (routes b and c; the dotted line is an interior crisis line) 
and boundary (routes d and e) crises may occur. Here the 
torus and the doubled torus are denoted by T and 2T, and 
the solid line represents a torus doubling bifurcation curve, 
(b) Smooth torus inside an absorbing area with boundary 
formed by segments of the critical curves Lk (k — 1, . . . , 5) 
for a = 3.38 and e = 0.584 7. (c) SNA filling the absorbing 
area for a = 3.38 and e = 0.584 75. For other details, see the 
text. 

determined by the critical curves. 

Figure 0(a) shows a phase diagram in the a — e plane. 
Each phase is characterized by the Lyapunov exponent 
a x in the x-direction as well as the phase sensitivity ex- 
ponent 5. The exponent S measures the sensitivity with 
respect to the phase of the quasiperiodic forcing and was 
introduced in || to characterize the strangeness of an at- 
tractor. A smooth torus that has a negative Lyapunov 
exponent without phase sensitivity (6 — 0) exists in the 
region denoted by T and shown in light gray. Upon 
crossing the solid line, the smooth torus becomes un- 
stable and bifurcates to a smooth doubled torus in the 
region denoted by 2T. Chaotic attractors with positive 
Lyapunov exponents exist in the region shown in black. 
Between these regular and chaotic regions, SNAs that 
have negative Lyapunov exponents with high phase sen- 
sitivity ((5 > 0) exist in the regions shown in gray and 



dark gray. Consistent with their positive phase sensitiv- 
ity exponent S, these SNAs are observed to have fractal 
structure Q . Here we restrict our considerations only to 
the intermittent SNAs that exist in the thin gray region 
[e.g., see a magnified part in Fig. 0(a)]. (In the dark- 
gray region, nonintermittent SNAs, born through other 
mechanisms, such as gradual fractalization JTcfl and torus 
collision [0], exist.) Note that this phase diagram is typ- 
ical for quasiperiodically forced period-doubling systems 
||, [13], g |T5|, [L6|, |lj|, and its main interesting feature 
is the existence of the "tongue" of quasiperiodic motion 
that penetrates into the chaotic region and separates it 
into upper and lower parts. We also note that this tongue 
lies near the terminal point (denoted by the cross) of 
the torus doubling bifurcation curve. When crossing the 
upper boundary of the tongue, a smooth torus is trans- 
formed into an intermittent SNA that exists in the thin 
gray region. Hereafter this intermittent route to SNAs 
will be referred to as the route a [see Fig. 0(a)]. 

As an example, we consider the case a = 3.38. Fig- 
ure 0(b) shows a smooth torus with a x — —0.059 for 
e = 0.584 7 inside an absorbing area whose boundary 
is formed by segments of the critical curves Lfc (k = 
1,...,5) (the dots indicate where these segments con- 
nect). We also note that the smooth unstable torus x = 
and its first preimage x = 1 form the basin boundary 
of the smooth torus in the 6 — x plane. However, as e 
passes a threshold value e* (= 0.584 726 781), a transition 
to an intermittent SNA occurs. As shown in Fig. 0(c) 
for e — 0.584 75, the newly-born intermittent SNA with 
<j x = —0.012 and 8 = 19.5 appears to fill the absorbing 
area, and its typical trajectory spends most of its time 
near the former torus with sporadic large bursts away 
from it. This intermittent transition may be expected to 
have occurred through collision of the smooth attracting 
torus with an unstable orbit. However, the smooth un- 
stable torus x — cannot interact with the smooth stable 
torus, because it lies outside the absorbing area. Hence 
we search inside the absorbing area for an unstable orbit 
that might collide with the smooth stable torus. 

Using RAs we find a new type of ring-shaped unsta- 
ble set that causes the intermittent transition through 
collision with the smooth torus. When passing the 
dashed curve in Fig. 0(a), such a ring-shaped unstable 
set appears via a phase-dependent saddle-node bifurca- 
tion. This bifurcation has no counterpart in the unforced 
case. (The dashed line is numerically obtained for the 
sufficiently large level k = 10 of the RAs.) For each 
RA of level fc, a periodically forced logistic map with the 
rational driving frequency u>k has a periodic or chaotic 
attractor that depends on the initial phase 6*o of the ex- 
ternal force. Then the union of all attractors for dif- 
ferent 6o gives the fcth approximation to the attractor in 
the quasiperiodically forced system. As an example, con- 
sider the RA of low level k = 6. As shown in Fig. 0(a) 
for a — 3.246 and e = 0.446, the RA to the smooth 
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FIG. 2: Ring-shaped unstable sets in the RA of level 6 for (a) 
a = 3.246 and e = 0.446, (b) a = 3.26 and e = 0.46, and (c) 
a = 3.326 and e = 0.526. (d) Ring-shaped unstable set in the 
RA of level 8 for a = 3.326 and e = 0.526. These ring-shaped 
unstable sets exist inside the absorbing area with boundary 
formed by portions of the critical curves Lk (k = 1, . . . , 4) 



torus, consisting of stable orbits with period F§ (= 8), 
exists inside an absorbing area bounded by segments of 
the critical curves Lk (k = 1 , . . . , 4) . Note also that 
a ring-shaped unstable set, born via a phase-dependent 
saddle-node bifurcation and composed of 8 small rings, 
lies inside the absorbing area. At first, each ring con- 
sists of the stable (shown in black) and unstable (shown 
in gray) orbits with the forcing period F 6 (= 8) [see the 
inset in Fig. |2|(a)]. However, as the parameters increase, 
a chaotic attractor appears through period doublings of 
the stable periodic orbit, and then it disappears through 
collision with the unstable Fg-periodic orbit (see Fig. ||(b) 
for a = 3.26 and e = 0.46). With further increase in the 
parameters, both the size and shape of the rings change, 
and for sufficiently large parameters, each ring consists 
of a large unstable part (shown in gray) and a small at- 
tracting part (shown in black), as shown in Fig. ||(c) for 
a = 3.326 and e = 0.526. [Note that the unstable part 
(plotted in the gray) of each ring consists of only un- 
stable orbits whose period is the same as that of the 
forcing, i.e., Fk-] For the same parameter values as in 
Fig. |2|(c), we increase the level of the RA to k = 8. Then 
the number of rings (= 336) increases significantly, and 
the unstable part [shown in gray and consisting of un- 
stable orbits with period Fg (= 21)] becomes dominant, 
because the attached attracting part (shown in black) be- 
comes negligibly small [see Fig. ||(d)]. In this way, with 
increasing the level k the ring-shaped unstable set con- 
sists of a larger number of rings with smaller attracting 
part. Hence, in the quasiperiodic limit these ring-shaped 
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FIG. 3: (a) and (b) RA to the smooth torus (shown in black) 
and the ring-shaped unstable set (shown in gray) in the RA 
of level 8 (F 8 = 21) for a = 3.38 and e = 0.586. (c) and (d) 
8th RA to the SNA for a = 3.38 and e = 0.586 4. For more 
details, see the text. 

unstable sets seem to form a complicated unstable set 
composed of only unstable orbits. 

In what follows we use RAs to explain the mechanism 
for the intermittent transition occurring in Figs. 0(b) and 
|l](c) for a = 3.38. Figures 0(a) and||(b) show that, inside 
the absorbing area, the ring-shaped unstable set (shown 
in gray) lies very close to the smooth torus (shown in 
black) for s = 0.586 in the RA of level k = 8. As 
e passes a threshold value £s (= 0.586 366), a phase- 
dependent saddle-node bifurcation occurs via a collision 
between the smooth torus and the ring-shaped unstable 
set, and then "gaps," where no orbits with period Fg, 
(= 21) exist, are formed. A magnified gap is shown in 
Fig. |(c) for e = 0.584 6. Note that this gap is filled 
by intermittent chaotic attractors together with orbits 
with period higher than Fg embedded in very small win- 
dows. As shown in Fig. ||(d), the RA to the whole at- 
tractor consists of the union of the periodic component 
and the intermittent chaotic component, where the lat- 
ter occupies the 21 gaps in and is vertically bounded 
by segments of the critical curves Lk (k = 1,...,5). 
However, the periodic component is dominant: the av- 
erage Lyapunov exponent (< a x > = —0.09) is negative, 
where < • • • > denotes the average over the whole 8. We 
note that Fig. ||(d) seems to be similar to Fig. |l|(c), al- 
though the level k = 8 is low. Increasing the level up 
to k = 15 we find that the threshold value ek at which 
the phase-dependent saddle-node bifurcation occurs con- 
verges to the quasiperiodic limit e* (= 0.584 726 781) in 
an algebraic way, |Aefe| ~ F k ~ a , where Aek — £k — £* 
and a ~ 2.2 ]2l| ]. Furthermore, in the quasiperiodic 
limit k —> oo the RA to the attractor has a dense set 
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of gaps that are rilled by intermittent chaotic attractors 
and bounded by portions of the critical curves. Thus, an 
intermittent SNA, containing the ring-shaped unstable 
set and filling the absorbing area, appears, as shown in 
Fig. |l|(c). In addition, we find that as e passes another 
threshold value e c (= 0.584 8), this SNA is transformed 
into a chaotic attractor with a positive Lyapunov expo- 
nent. Using the RA, this transition to chaos may also be 
explained. For each RA, its angle averaged Lyapunov ex- 
ponent < a x > is given by the sum of the "weighted" Lya- 
punov exponents of the periodic and chaotic components 
of the RA, A p and A c , (i.e. < a x >= A p + A c ), where 
A p(c) = Mp( c ) < a x >p( c ), and M p(c ) and < a > p(c) are 
the Lebesgue measure in 9 and average Lyapunov expo- 
nent of the periodic (chaotic) component, respectively. 
After passing a threshold value where the magnitude of 
A p and A c are balanced, the chaotic component becomes 
dominant, and hence a chaotic attractor appears. (More 
details will be given in Ref. |2l|.) 

To sum up, using RAs we have found the mechanism 
for the intermittent route to SNAs in the quasiperiodi- 
cally forced logistic map. When a smooth torus makes 
a collision with a new type of ring-shaped unstable set, 
a transition to an intermittent SNA, bounded by seg- 
ments of critical curves, occurs via a phase-dependent 
saddle-node bifurcation. Other typical transitions such 
as interior [routes b and c in Fig. |l|(a)] and boundary 
[routes d and e in Fig. 0(a)] crises may also occur near the 
main tongue through interaction with the ring-shaped 
unstable set Furthermore, as e decreases toward 

zero, similar tongues appear successively near the termi- 
nal points of the higher-order torus-doubling bifurcations 
p2j, and band- merging crises may also occur through 
collision with the ring-shaped unstable set pl| . Conse- 
quently ring-shaped unstable sets play a central role for 
dynamical transitions occurring near the tongues. Fi- 
nally, we note that these kinds of dynamical transitions 
seem to be "universal," in that we observe that they 
occur in typical quasiperiodically forced period-doubling 
systems of different nature, such as the quasiperiodically 
forced Henon map, ring map, and pendulum |21|] . 
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